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i=n 

2 (x{ cos a + yi sin a) 2 = Minimum. 
Equating to zero the derivative as to a, 

S {y? ~ x ?) sm a cos a = 2 ^i2/i(sin 2 ce — cos 2 a). 
Let 
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The origin shall be 0, the point (1, 0) shall be 0'. We construct a unit circle 
about 0' as center. Let OPi meet this circle again in Q t . The coordinates of 
Qi in the system w/iere 0' is origin will be (cos 20»-, sin 20*). Let Q be the center 
of gravity of the points Qi. The line O'Q shall meet the circle in R, that is to 
say, R is the intersection which lies on the same side of the diameter 00' as does 
Q. Then OR is the required line. See Fig. 2. 



COMPUTATION FORMULA FOR THE PROBABILITY OF AN EVENT 
HAPPENING AT LEAST C TIMES IN N TRIALS. 

By B. C. MOLINA, New York City. 

Many problems in biometry, radioactivity, etc., require for their quantitative 
solution a convenient formula for computing the probability of an event happen- 
ing at least c times in n trials; the probability p of its happening in one trial 
being known. Several relatively simple formulas are given by Poisson in the 
third chapter of his Recherches sur la ProbabiliU des Jugements, availing himself 
of a method developed by Laplace in the TMorie Analytique des Probability. 
But these formulas explicitly exclude the case where p is very small and 
c = np + r, r not being small compared with np. For this range of values a 
problem connected with my engineering work forced me early in 1908 to develop 
the formulas* (5) and (6) given below. 

Let P = the required probability, a = np, s = (n — c)/(c + 1), 

2e(l — x)x 2 + (n — s)sx A 



F(x) = 



[c — in — s)x] i 



* These formulas have been used for the construction of curves which are in constant use 
in the Engineering Department of the American Telephone and Telegraph Company. They are 
published in the hope that others will find them helpful. 
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and take as a starting point Poisson's exact equation 

Ydy y n ~° 1 - v 

f = -= , where Y = . )n+1 and <r = — — . 

J Ydy K ^ y ' V 



Transforming the numerator by the change of variable y = (l/x) — 1 and 
substituting for the integral in the denominator its well-known value, gives 

ni r 

(1) P = -. ^-77 n xT\l - x) n ~ c dx, 

(c — l)i(ii — c)\J 

or, as 

c = np + r > in — s)p, 

we may write 

P = — C — \x°(l — x) n - c ~ s \ — — X ' >S+1 dx 

r (c - 1)1 (» - c)l Jo dx iiru X) } [c - (n- s)x] aX - 

By partial integration 

, 9 v p nl f p°{l - p) n ~° +1 _ p' gHq - xT2 1 

W ^ _ (c-l)!(n-c)!l[c-(n- 5 )p] lra 5;s J [e-(n - s)xf dX | " 

Also, by partial integration 

r x° +1 (i - x) n ~° __ r e_ x 2 {i - x) s+1 

J [c-(n- s)xf dX ~ X dx {X ° {1 xr ' } [0 - (n - s)xf dx 
(3) 

p c+2n _ v \n-c+l ,'P 

[c — (ra — s)p] 3 Jo 

Now, since. P(a;) remains finite, does not change sign and is a maximum for 
a; = p as a; varies from to p, provided p < 2/3, 

a^-^l - x) n ~ c F(x)dx = \F(p) xT\l - x) n ~ c dx, < X < 1. 
Jo 

By (2), (3), (4) and the expressions for s and F(p) we obtain 



P nl en y~f (l-y)(c+D 1 



c! (n 



(5) 



c(c + 1 — oc — p) 2 
where 



M 



\cp 



(a + p)(a- cp) {2(c + 1) 2 (1 - p) + (a + p) (a - op) } o 2 {2(c + l) 2 + c?\ 
V= „2/„ 11 „m < 



c 2 (c +l-a-p) i ^ c?(c+ 1 - a - p) 4 
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For p not greater than 1/30 and P not greater than .01 I found that the values 
of c and a (= np) are such that <p is negligible. For this range of values, therefore, 
the factor 1/(1 — \<p) in the right-hand member of (5) was ignored and a con- 
venient formula for computing P thereby obtained. After plotting several 
curves it became evident that for a given value of P the corresponding value of c 
did not change appreciably with changes in the values of p and n, provided the 
product a (= np) was kept constant. The following formulas were therefore 
used. 

In (5) let n and p approach the limits oo and respectively while a remains 
constant. Then in the limit 



1 1 I ^£2 1 I c+1 1 
-X^jj c! ]\e+ 1-aJ 



a' 



1-XlH c! JU+l-aJl c(c+l-a) 2 

(6) 

a 2 {2(c+l) 2 + a 2 } 

*'" c^c+l-a) 4 ' 

Let P , Pi be the values of P obtained from (6) by taking X equal to and 1 
respectively. Then 

(7) Po < P < Pi. 

Therefore (Pi + Po)/2 may be taken as an approximate value of P provided 
(P L — P )/2 is small compared with P . This will be the case if P is not greater 
than .01 and a/c small. As a/c increases the discrepancy between P and 
(P ± + Po)/2 increases, but not indefinitely. This may be proved as follows. 

For a given value of P we know (Bernoulli-Laplace Theorem) that a/c ap- 
proaches the limit 1 as a and c approach oo . Also, for a and c very large, one 
of Poisson's formulas is a very close approximation and reduces to 

2 r a 

(8) P = ^= e ~^ 

where k, the lower limit of the integral, is given by the equation 
k 2 + (c - 1 - a) 



(c-1) 



= - log 



a 1 . j c—l—a 
— --log 1 ;:=T - 



(c - 1 - a) (c - 1 - ay (c-l- a y 
~ (c-1) + 2(c - l) 2 + 3<c - l) 3 

or, for c very large and a/c = 1 



(9) «« - 
but for a/c = 1 

(10) 4> = 



, _ (c - 1 - a) 2 _ (c + 1 - a) 2 



2(c - 1) 2a 

but for a\c = 1 

a 2 {2(c + l) 2 + a 2 } 3a 2 



c»(o + 1 - a) 4 (c+1 - a) 4 (2k 2 ) 2 
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Therefore, by (10) and any of the well-known tables giving the values of P and k 
which satisfy equation (8), we obtain the following figures showing the limits 
(as c = oo and a/c = 1) of the possible error incurred by taking P = (Pi + Po)/2. 



p 


K 2 


* 


(Pi - Po)/2P» 


.01 


2.72 


.101 


.056 


.001 


4.80 


.032 


.017 


,0001 


6.92 


.016 


.008 



The short table which follows gives some of the pairs of values of c and a satis- 
fying the three equations P = .01, P = .001 and P = .0001 for the limiting 
case of p = 0, n = °o . 





P = . 


,0001 


P = , 


,001 


P = 


.01 


c 


a 


a/c 


a 


a/c 


a 


a/c 


1 


.0001 


.0001 


.0010 


.001 


.010 


.010 


2 


.0142 


.0071 


.0454 


.023 


.149 


.075 


3 


.0862 


.0289 


.191 


.064 


.436 


.145 


4 


.232 


.0580 


.429 


.107 


.823 


.206 


5 


.444 


.0888 


.739 


.148 


1.28 


.256 


6 


.714 


.119 


1.11 


.185 


1.79 


.298 


7 


1.03 


.147 


1.52 


.217 


2.33 


.333 


8 


1.39 


.174 


1.97 


.246 


2.91 


.364 


9 


1.78 


.198 


2.45 


.272 


3.51 


.390 


10 


2.20 


.220 


2.96 


.296 


4.13 


.413 


12 


3.11 


.259 


4.04 


.337 


5.43 


.452 


14 


4.11 


.293 


5.20 


.371 


6.78 


.484 


16 


5.17 


.323 


6.41 


.401 


8.18 


.512 


18 


6.28 


.349 


7.66 


.426 


9.62 


.534 


20 


7.44 


.372 


8.96 


.448 


11.1 


.555 


30 


13.7 


.457 


15.9 


.530 


18.7 


.623 


40 


20.6 


.515 


23.3 


.582 


26.8 


.670 


60 


35.4 


.590 


38.0 


.648 


43.5 


.725 


80 


50.9 


.636 


55.2 


.690 


60.7 


.759 


100 


67.0 


.670 


71.9 


.719 


78.2 


.782 


140 


100.2 


.716 


106.3 


.760 


114.0 


.814 


180 


134.3 


.746 


141.4 


.786 


150.3 


.835 


200 


151.6 


.758 


159 


.795 


169 


.845 


300 


240 


.800 


249 


.830 


261 


.870 


400 


330 


.825 


341 


.852 


355 


.887 


500 


421 


.842 


434 


.868 


450 


.900 


600 


513 


.855 


527 


.878 


545 


.908 


700 


606 


.866 


621 


.887 


640 


.914 


800 


699 


.874 


715 


.894 


736 


.920 


900 


793 


.881 


810 


.900 


832 


.924 


1000 


887 


.887 


905 


.905 


928 


.928 



BOOK REVIEWS. 

W. H. BUSSEY, Chairman of the Committee. 

Academic Algebra. By George Wentworth and David Eugene Smith. 

Ginn & Co., Boston, 1913. iv + 442 pages. $1.20. 

This book is the second of the Wentworth-Smith series. In it the authors 
have attempted to provide a high school course which shall cover the topics 



